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Abstract 

We have measured the cross section a{e^e~^TT~^TT~^) with the KLOE detector 
at DA<1>NE, at an energy W = = 1.02 GeV. Prom the dependence of the 

cross section on m(7r+7r~) = ^Jw^'^-2WE^, where E-y is the energy of the photon 
radiated from the initial state, we extract (T(e"'"e~^7r+7r^) for the mass range 0.35 < 
m2(7r+7r-) < 0.95 GeV^. From our result we extract the pion form factor and the 
hadronic contribution to the muon anomaly, a^. 



1 Hadronic Cross Section at DA$NE 

1.1 Motivation 

The recent precision measurement of tlie muon anomaly at tlie Brookhaven 
National Laboratory [1] has led to renewed interest in accurate measure- 
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ments of the cross section for e^e~ annihilation into hadrons. Contributions 
to the photon spectral functions due to quark loops, are not calculable for 
low hadronic mass states because of the failure of perturbative QCD in such 
conditions. A very clean way out has been known for a long time [2]. The 
imaginary part of the hadronic piece of the spectral function is connected by 
unitarity to the cross section for e"''e~ ^hadrons. A dispersion relation can 
thus be derived, giving the contribution to cis ciri integral over the hadronic 
cross section multiplied by an appropriate kernel. An example of a first com- 
plete estimate of the correction was given in 1985, 6 a^^'^ =7 07 (6) (17) x 10~^° 
[3]. The process e~^e^^7i^7i^ contributes ~500 out of the ~700 value above. 
The cross section for e^e~-^7r+7r~ becomes negligible above 1 GeV. 

The most recent measurements of ale^e^^n^n") in the energy interval 610< 
m{n^7r~) <961 MeV come from CMD-2 at Novosibirsk. They claim a system- 
atic error of 0.6% with a statistical error of ~0.7% [4]. Their data have been 
used most recently together with r and e+e~ data up to 3 GeV, in an attempt 
to produce a firm prediction for comparison with the BNL result [5]. There is 
unfortunately a rather strong disagreement between the 6a^^'^ value obtained 
using e"'"e~— >7r+7r~ data and r decay data, after isospin corrections. Finally the 
e"'"e~— ^TT+TT" based result disagrees by ~3 cr with the BNL [1] measurements. 
There are thus many reasons for new measurements of the e~^e~ annihilation 
cross section into two pions. 



1.2 Initial State Radiation 

Initial state radiation, ISR, is a convenient mechanism which allows studying 
e"'"e~— s> hadrons, over the entire energy range from 2mTr to W, the center of 
mass energy of the collision. In the case of interest it is potentially vitiated by 
the possibility of final state radiation. For a photon radiated before annihila- 



tion of the e~^e~ pair, the vr'^vr" system energy is m(7r+7r~) = — 2WE^, 

thus one measures the coupling to pions of an off-mass shell photon of mass 
m(7r+7r~). For a photon radiated by the final state pions, the virtual photon 
coupling to the tt+tt" pair has a mass W. By just counting powers of a, the 
relative probability of ISR and FSR are equal. This requires very careful esti- 
mates of the two processes in order to be able to use the reaction e~^e~^7i^7i~'y 
to extract a{e~^e^^7r~^'iT~). The Karlsruhe group has given much attention to 
this problem and has also developed the EVA and PHOKHARA Monte Carlo 
programs [6] - [10] which are fundamental to our analysis. In the following we 
will refer to the invariant mass squared of the tt+tt^ system as s^. In general 
the 7r"'~7r~7 and tt'^tt^ cross section are related through: 




da"(7r'^7r 7) 



s. 



•n 



dsTT d cos 9 



= (7(71 71 
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Integrating the cross section above, for < 6*^ < 6* and 180° — 6 < 6^ < 180 
we get 

Sn — -\ = (T(vr+7r , s^) X H{s^, 9) (1) 



Eq. 1 defines tlie radiator function H{s.^, 9). In the following we will drop the 
variable 9 which is a constant in the present work. The radiator function H 
used in our analysis is obtained from the PHOKHARA Monte Carlo program. 
Our present analysis is based on the observation of ref.[6], that for small 
polar angle of the radiated photon, the ISR process vastly dominates over the 
FSR process. At lowest order the tx'^tx^^ cross section diverges as 1/ sin^ ^, 
just like (T(e"'"e~— +77). We limit ourselves in the following to studying the 
reaction e~ ^ti^ t[~ '-y with 9^ < 15° or 9^ > 165°. For small m(7r"'"7r~), the 
di-pion system recoiling against a small angle photon will results in one or both 
pions being lost also at small angle. We are therefore limited to measuring 
a^ir^TT^) for ■m{7r~^7T~) >550 MeV. We will be able to investigate the cross 
section near threshold, as soon as next to leading order calculation for FSR 
become available. This is of great importance, since there are no recent, good 
measurements of (j(-7r"'"7r~) at low mass, which weigh strongly in the estimate 
of 5a^<^. 



2 The KLOE detector 



The KLOE detector [11] consists of only two elements: a large, precision drift 
chamber and an electromagnetic calorimeter. The detector design was opti- 
mized for CP study in neutral kaon decays, for K mesons produced in the 
decay of mesons almost at rest, P(^~12.5 MeV/c. The measurements re- 
ported in the following rely heavily on the drift chamber. Its large dimensions, 
4 m diameter, use of He and a magnetic field of 0.52 T allow measuring mo- 
menta and direction of charged particles with very high accuracy. Because 
of our choice of accepting only events with |cos6'^| > cos 15°, the photon 
does not reach the calorimeter. The two relevant parameters, and 9^ are 
obtained from the reconstructed tracks of events with two opposite charge 
particles forming a vertex near the e'^e~ interaction point. The quality of the 
detector is fully described by the mass squared resolution and the resolution 
in measuring the production angle of single tracks or the angle of the vector 
sum of the two track momenta. We find cr(m(7r'^7r~))/m~0.2% for m=500 
MeV, without much change over the range studied and (t(^^)~10 mrad, with 
9^ = cos~^(p_,__ ■ z) and p+_ = — (p"*" + P~)/|(p''' + P~)|- The momentum 
resolution is (T{p±)/p± < 0.4% for large angle tracks. The calorimeter still 
plays an important role in the measurement. Its exceptional timing resolution. 
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a{t) = 57 /VE GeV © 50 ps and high segmentation provide particle identifi- 
cation allowing us to distinguish pion from electrons and muons, as discussed 
later. It is also used in the luminosity measurement, 

The KLOE detector is operated at DA$NE, the Frascati 0-factory, running 
aX W = = 1019.4 MeV. Beams cross in DA$NE at an angle of 25 mrad, 
resulting in a momentum of the mesons of ~12.5 MeV/c along the x-axis. We 
use a coordinate system where the a;-axis points to the center of the colliders, 
the z-axis bisects the two beam lines and the y-axis is vertical. The precise 
values of W and is measured on a run by run basis using Bhabha events, to 
accuracies significantly better than 100 keV. The rms machine energy spread 
is ~350 keV, measured with KsK^ pairs in KLOE. 



3 Event Selection 



The excellent KLOE resolution allows us to efficiently selected 7r^7r^7 events, 
without photon detection, because of their clean signature: two tracks from 
the interaction point. We assume that the event is of the type e"''e~^x"'"x~7, 
with m{x^) = m{x~) and = m^. The undetected photon polar angle, 6^ 
is the angle with the ^-axis of the vector -p+-p-. 



3.1 Fiducial Volume 



The fiducial volume for 7r"'"7r~7 events has been chosen in order to maximize 
the acceptance for 7r'^7r~7 events due to ISR. The main background processes 
are e"'"e~— s-7r"'"7r~7r° and radiative Bhabha scattering. At low mass the final 
state /i"'"/i~7 becomes important but requiring small 6^ removes most of this 
final state. We rely otherwise on kinematics and particle identification, ID, in 
order to reject the mentioned backgrounds to a negligible level. From studies 
with EVA, confirmed with PHOKHARA, a good compromise for the fiducial 
volume is | cos6'^| > cos 15° and 50° < 6'_|__ < 130 for the polar angle of both 
tracks at the interaction point. These requirements ensure abundant statistics 
for the measurement and contamination of FSR hard radiation well below 
the 1% level. The efficiency for the events falling inside the fiducial cuts is 
contained in the value of the radiator, obtained from PHOKHARA. We note 
that no detected photon is required, nor are events rejected if photons are 
present. 
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3.2 Event Selection 
3.2.1 Event origin 

The two tracks must form a vertex close to the interaction point. We require 



that the vertex coordinates a:y,Xy,X2 satisfy + < 8 cm and \zy\ < 7 
cm. We also require Px_ > 160 or > 90 MeV/c to remove spiraling tracks. 90 
to 95% of the 7r"'"7r~7 events survive these requirements. This result is obtained 
from Monte Carlo studies and examination of 0^7r"'"7r~ decays. 



3.2.2 Pion identification 

A fraction of radiative Bhabha events, e^e~7, survives kinematical selection, 
see sec. 3.2.3 below, contributing a non-negligible background. In order to 
reject these events, we have developed a particle ID method, using approximate 
likelihood estimators which are effective in distinguishing pion and electrons. 
The likelihood is based on time-of-flight vs momentum comparison and on 
the shape and magnitude of the energy deposits in the calorimeter by the 
charged particles. Correlations, which are significant, between used variables 
are ignored. Two likelihood functions have been built using actual data: one, 
to test the pion hypothesis^ (using -7r"''7r~7r*^ events) and the other, to 
test the electron hypothesis (using Bhabha events). We require that candidate 
pion tracks have a greater likelihood of being pions than of being electrons. We 
use as discriminating variable ( = log(L^/Le) and we accept as pions particles 
with C > 0, see fig. 1. The efficiency for pions is ~98% and 3% of the electrons 
survive rejection. We require that only one of the two particles satisfy the pion 
ID requirement (OR) achieving an efficiency of 100% with some very small 
residual background. Requiring that both tracks are pions (AND) leads to a 
5% loss. We use the AND for evaluating backgrounds, see later. 



3.2.3 Kinematic Cuts 

For an x~^x~'~f final state Lorentz invariance allows us to compute the mass of 
X, which we call mtrk from 



^ In fact, in the case of the pion likelihood, two different functions are built, one 
for 7r+ and one for vr", since the energy deposition of the two particles is different, 
in particular at these energies. This is not necessary for electrons. 
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Fig. 1. Distribution of the particle ID variable C for pions, right peak, and electrons, 
left. See text. 

For e+e~ x~^x~'-f, the value of m-trk peaks at ttItt, m^, rrie thus allowing 
selection of signal events. The density distribution of the two track events in 
the [Stt, mtrk] plane is very effective for separating signal from background. 
The distribution is shown in fig. 2 together with the selection cut defined by: 
(mtrk > 120) © (mtrk < 250 - 105^1 - (mtrk < 220), all units in 

MeV. 7i~^7r~7T^ events populate a curve in the plane. The signal events above 
the region mtrk > 140 MeV are due to additional ISR or FSR radiation. 
The cut on mtrk rejects a fraction of events with additional radiation. The 
implications are discussed in section 4.4. 



4 Event Analysis 



Data were taken in July through December 2001, for an integrated luminosity 
of L=140.7 pb~^. After fiducial volume and selection cuts we find ~1.5 x 10^ 
events. Fig. 3 gives the vr"*"?!" mass squared distribution of the 7r'^7r~7 events 
in bins of 0.01 GeV~^ in s^. The p peak and the p-u interference structure 
are clearly visible, demonstrating the excellent mass resolution of the KLOE 
detector. To obtain the cross section, we subtract the residual background 
from this spectrum and divide by the selection efficiency and the integrated 
luminosity: 

Acr^+^-^ _ AA^+^-^ _ AApbs — AAbckgnd ^ 1 
As^ L As^ L As^ ecuts 
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4-1 Background 



The [st,, rritrk] plane population of signal and background events is illustrated 
in fig. 2. The amount of background in the signal region is obtained by fitting 
the rritrk spectrum of the selected events, in slices of s^, with the mtrk spectra 
for signal, e+e~7 and fi'^fi~'-f events obtained from Monte Carlo simulation. 
An example of such a fit is shown in fig. 4, left. Background from 7r"'"-7r~7r° 
events appears at higher mtrk values and the missing mass, m^j^^ = {p^ — 
— peaks at m(7r°). The tt^tt^tt^ background count is obtained by 

fitting the mmiss distribution with the shapes obtained from the Monte Carlo 
simulation. An example of such a fit is shown in fig. 4, right. The result of 
the fits is shown in fig. 5 which shows the relative contribution of the different 
backgrounds to the observed signal. The shape of the background distribution 
is well reproduced by the Monte Carlo simulation ensuring that systematic 
uncertainties are smaller than the fit errors. Other possible contributions, such 
as the reaction e"'"e~— *-e"'"e~7r"'"7r~ with the electrons going down the beam pipe, 
have been studied and shown to be negligible given our event selection criteria. 



4.2 Luminosity Measurement 



The integrated luminosity is measured using large angle Bhabha (LAB) events, 
with the KLOE detector itself. The effective Bhabha cross section at large 
angles (55° < < 125°) is about 430 nb. This cross section is large enough 
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Fig. 2. Kinematic separation, in the [rritrk, ■Stt] plane, of signal and backgrounds after 
the C variable cut. 
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Fig. 3. Events per 0.01 GeV after fiducial volume and selection cuts. 




Fig. 4. Muon mass fit, left and vr*^ mass fit, right. These fits are used to estimate 
the and 'k^'k~'k^ backgrounds (shaded areas) to the 7r"'"7r~7 channel. 

so that the statistical error on the luminosity measurement is negligible. The 
number of LAB candidates, N^aB) is counted and normalized to the effective 
Bhabha cross section, obtained by Monte Carlo: 



/ 



= MC (n ^ ■ (1 - ^Bkg) (4) 



The precision of this measurement depends on the correct inclusion of higher 
order terms in computing the Bhabha cross section. We use two independent 
Bhabha event generators: BHAGENF ([14], [15]) and BABAYAGA ([16]). For 
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Fig. 5. Left: backgrounds fractions requiring a single identified pion (OR) and right: 
two identified pions (AND) which removes completely the residual Bhabhas. 

each generator a systematic error of 0.5% is quoted by the authors. The two 
result for the luminosity agree to better than 0.2%. 

LAB events are selected with cuts on variables which are well simulated by the 
KLOE Monte Carlo simulation. The electron and positron polar angle cuts, 
55° < 6^^^ < 125°, is based on the calorimeter clusters, while the energy cuts, 
E'^ _ > 400 MeV, is based on drift chamber information. The background from 
/i+yU~(7), 7r"'"7r~(7) and 7r"'"7r~7r° events is well below 1% and is subtracted. All 
selection efficiencies (trigger, EmC cluster, DC tracking) are > 99% and are 
well reproduced by the detector simulation program. We also obtain excellent 
agreement between the experimental distributions (6'_|_^_, -£'+,-) and those ob- 
tained from Monte Carlo simulation, see fig. 6. Finally, corrections are applied 
on a run-by-run basis for fiuctuations in the center-of-mass energy of the ma- 
chine and in the detector calibrations. The experimental uncertainty in the 
acceptance due to all these effects is 0.4%. We assign a total systematic error 
for the luminosity of SL = 0.5%th © 0.4%cxp- An independent check of the lu- 
minosity measurement is obtained using e"'"e~— ^77 events. We find agreement 
to within 0.2%. 



4-3 Selection Efficiencies 



4-3.1 Trigger, tracking and vertexing 

These efficiencies have been obtained from control samples of 7r"'"7r~7r'^, vt+tt" 
and 7r"'"7r~7 events selected without requiring the vertex. The efficiencies have 
been parameterized as a function of single track parameters (charge, momen- 
tum and angle) combined using Monte Carlo simulation. The largest source 
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Fig. 6. Comparison data-Monte Carlo of the 0+,- (left) and -£'+,- (right) distribu- 
tions for the Bhabha events selected at large angle as described in the text. 



of systematic error in this procedure arises from the difference between the 
kinematics of the control samples and the signal. 



4-3.2 Pion identification 

. The efficiency for the C, cut has been evaluated from data by accepting one 
track and studying the distribution of the other one. The systematic error 
introduced by the C, cut is negligible, since we use the OR, i.e. at least on 
particle is a pion, for which the efficiency is ~100% for 7r'^7r~7 events. 



4-. 3. 3 Particle mass 

The efficiency for the mtrk cut is obtained as a by-product of the residual 
background evaluation; the result of the fit provides the efficiency in each s.,, 
bin. However, this efficiency depends upon the way the Monte Carlo simula- 
tion treats multi-photon processes. The mtrk efficiency has been obtained with 
our reference Monte Carlo simulation, which uses the PHOKHARA version 
1.0 generator. PHOKHARA simulates the full NLO QED correction for e"^e~ 
annihilation with the emission of additional photons from the initial state. In 
order to estimate the systematic effect due to the generator we have repeated 
the exercise using a different generator, BABAYAGA ver.3.5, the same genera- 
tor used in the luminosity measurement. In this generator, ISR is treated using 
the parton shower approach. The resulting value for the mt^k efficiency differs 
from that evaluated with PHOKHARA by 0.2%. An additional contribution, 
which is treated separately, is the emission of both an ISR and FSR photon. 
This is discussed in the next section. The contribution of the systematic error 
for each efficiency to the total error is shown in table 1. 
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Acceptance 


0.3% 


Trigger 


0.6% 


Tracking 


0.3% 


Vertex 


1.0 % 


Likelihood 


0.1 % 


Track Mass 


0.2 % 


Background subtraction 


0.5 % 


Unfolding 


0.6 % 


Total exp systematics 


1.4 % 


luminosity 


0.6 70 


Vacuum Polarization 


0.2 % 


Total theor. systematics 


0.7% 


FSR resummation 


2.0 % 



Table 1 

List of systematic uncertainties from the three sources: experimental, theoretical 
and ignoring FSR. 

4-4 Unfolding the mass resolution 

In order to obtain d(T,r7r7/ds7r as a function of the true s-,, value, we unfold 
the mass resolution from the measured s.^ distribution. The measured value of 
■57r,obs is related to the true value by a resolution matrix G(s7r,truc — s^^obslsTr.truc)- 
This matrix has been generated by Monte Carlo simulation. The resolution 
matrix is nearly diagonal with 75% of the events in the diagonal, 11% in each 
of the first once-off-diagonal elements and only 3% in the remaining elements. 
The matrix is graphically illustrated in fig. 7, left. Unfolding is obtained by 
applying the inverse matrix to the data. The result is shown in fig. 7, 
right, which shows the difference between the reconstructed and the generated 
spectra with and without the application of the unfolding procedure. This 
procedure does not affect the value of a^, since the correlation is short range, 
^ 97% of the events being in the same or the adjacent bins. 



5 Results 

The cross section for e^e" 7r+7r~7, after applying the corrections de- 
scribed above, is shown in fig. 8. In order to get the e^e" tt+tt^ cross 
section, the radiator function H{st^^ 9^) is needed. The radiator is obtained 
from PHOKHARA, setting F^(s^) = 1 and switching off the vacuum polariza- 
tion of the intermediate photon in the generator. The radiator H{s.^, 9^ = 15°) 
is also shown in fig. 8, left. 
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Fig. 7. Left: smearing matrix representing the correlation between generated and 
reconstructed value; the high precision of the DC results in an almost diagonal 
matrix. Right: fractional difference between reconstructed and generated spectrum 
with and without resolution unfolding. 
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Fig. 8. Left. Cross section for e'^e tt^tt 7. The radiator H[st^, 6^), is also shown. 
Right. Bare cross section for e~^e~ — > 7r"*'7r~. 



5.1 FSR and Vacuum Polarization corrections 



Part of the multiphoton events are removed by the mt^k cut. The effect of this 
cut is evaluated using PHOKHARA 1.0 which, however, does not include soft 
radiation in the final state accompanying one hard photon from the initial 
state. An preliminary estimate done with PHOKHARA-11, shows that this 
effect is smaller than 2%. We do not correct here for FSR photons and we 
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include in the theoretical systematic error the full value (2%) due to this effect. 
In order to obtain the pion form factor, or the bare cross section, vacuum 
polarization effects must be subtracted. This can be done by correcting the 
cross section for the running of a, see refs. [17], [5], as follows: 

'^bare ^dressed 



\a{s) 



We have used an approximate, but quite adequate for the purpose estimate of 
a(s)[18]. Fig. 8, right, shows the resulting cross section. 



5.2 Results and comparisons with existing data 



We have used the bare cross section to evaluate of 6ajf' in the region covered 
by the CMD-2 experiment (0.37 < < 0.93). The resulting value (in 10"^'^ 
units) is : 



Sal"^ = 374.1 ± l.lstat ± 5.2syst ± 2.6theo lao 



FSR 



(6) 



to be compared with: 

5aJ'^'^ (CMD-2) = 368.1 ± 2.6stat ± 2.2syst+the 



(7) 



The statistical error is negligible. The systematic error is at the moment at 
the 1.4% value, but there is still room to improve it down to ~ 1%. Clearly the 
dominating error is the one we have conservatively assigned to the FSR radia- 
tion effect. This error will be substantially reduced once the new PHOKHARA 
version will be inserted into the KLOE MC. 

The central value of our result in the region 0.37 < s^^ < 0.93 is slightly larger 
than the one obtained by CMD-2. The discrepancy is, at the moment, not 
significant (0.5 a). This difference is, however, not equally distributed in s,r, 
as summarized below. ^ 

s^, GeV^ 5a^, KLOE 6a^, CMD-2 

0.37 - 0.60 256.2 ±4.1(t^-^ FSR) 249.7 ±2.2 
0.60 - 0.93 117.9±2.1(lg-^ FSR) 119.8 ±1.1 

Our data differ from the CMD-2 results mostly below the p-peak. However, 
for the mass squared range 0.6 < < 0.9, our data confirm the discrepancy 
between r data and e^e~-^7r+7r^ results, which is ~10-15%. 

^ This is our evaluation of aj^'^'^ (CMD-2), based on the values tabulated in ref. [4] 
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